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Abstract 

Small particles transported by a fluid medium do not necessarily have to follow the flow. We 
show that for a wide class of time-periodic incompressible flows inertial particles have a tendency 
to spontaneously align in one-dimensional dynamic coherent structures. This effect may take 
place for particles so small that often they would be expected to behave as passive tracers and 
be used in PIV measurement technique. We link the particle tendency to form one-dimensional 
structures to the nonlinear phenomenon of phase locking. We propose that this general mechanism 
is, in particular, responsible for the enigmatic formation of the 'particle accumulation structures' 
discovered experimentally in thermocapillary flows more than a decade ago and unexplained until 
now. 



Ordering and transport of small particles suspended in an incompressible fluid medium 
is both fundamentally interesting and important in an array of natural and technological 
processes. Their examples include manipulation and segregation of biological and synthetic 
microparticles [Ij, concentration and transport of pollutants, nutrients, and plankton by 
oceanic currents and of aerosols and droplets by atmospheric currents [i2j. 

When interparticle interactions are negligible, small particles are often expected to be- 
have as passive tracers when their Stokes number, defined as the ratio of the particle re- 
sponse time to the characteristic fiow time scale, is much less than unity. On the other 
hand, it has been well-known in the fiuid mechanical literature that small but finite sized 
particles with densities different from that of the surrounding fiuid may deviate from the 
fiuid trajectories. These deviations, accrued over time, lead to non-conservative trajectories 
in otherwise volume-preserving fiows. This fundamental effect may have profound conse- 
quences. In the past it has been numerously invoked to explain segregation and unmixing of 
heavy particles [3j . Recently a strong theoretical interest to this effect resulted in predicting 
the possibility of permanent suspension of particles in fiows under gravity [1], formation 
of three-dimensional particulate coherent structures in hurricanes [5j, trapping of aerosol 
particles in open chaotic fiows [6j and in a chaotic wake [7J. A systematic treatment of the 
coherent structures arising due to inertia effects along the lines of dynamical systems was 
undertaken recently [5j. 

In this Letter we report a new type of ordering of inertial particles that results in formation 
of one- dimensional dynamical particulate coherent structures. It is rather surprising that 
while this effect is generic for a class of widely encountered time-periodic fiows, it has 
apparently not been analyzed previously. It is even more surprising that the inertia-driven 
self-assembly of particles into one-dimensional continuous lines was actually observed in 
experiments on thermocapillary fiows more than a decade ago by Schwabe et al. [8j but 
despite extensive experimental studies P, [lO] has remained unexplained. Our results were 
first reported at pT] . 

The fiuid fiows we consider occur in the cylindrical geometry. They can be represented 
as superpositions of a steady toroidal vortex u^^\r^z) and an oscillatory wave traveling 
azimuthally u^^\(/) — Qt^r^z)^ see Fig. [!} Here (0, r, z) are the cylindrical coordinates, t 
is time, and Q is the wave angular velocity. Such fiows often emerge as a result of an 
instability, e.g. the convective instability in a rotating annulus that had been suggested as a 
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FIG. 1: Sketch of a time-periodic flow with an azimuthal travehng wave. Here the wave mode 
m — 2. Such flows often occur in the cylindrical geometry as a result of an instability [121 [131 [T5]. 
In the experiments of Schwabe et al. the flow was induced by the Marangoni effect [16j in a drop 
placed between two rods, with the top rod heated. 

model for atmospheric circulation [12j or the instability in a stationary open cylinder driven 
by rotating endwall [13j. The ubiquity and practical importance of such flows is the reason 
why a number of studies have been devoted to them, particularly with the focus on chaotic 
advection and mixing [14j. Here, however, our focus is different. 

In the experiments of Schwabe the instability was of thermocapillary origin [15]: it devel- 
ops when a drop of liquid is placed between two cylindrical rods, with the top rod heated. 
The flow, driven at the cylindrical liquid-gas interface by the Marangoni force [16j, under- 
goes a bifurcation from a steady axially-symmetric vortex to an oscillatory regime with a 
traveling wave when the temperature difference between the rods is large enough. That 
particles will self-assemble in this flow was discovered when the experimentalists admixed 
particles, having sizes of tens of microns, in order to study the flow. To much surprise, they 
discovered that under certain conditions such small particles defied the fate of passive trac- 
ers and aligned themselves in an ordered spiral structure (see Fig. [2](a) and Supplementary 
Movie 1 [17j). The spiral was closed, rather symmetric and rotated around the axis with 
no change of shape. It was dubbed TAS' for 'particle accumulation structure' [8j. Notably, 
the angular frequency of rotation of the spiral was found equal to that of the wave. In other 
experiments particulate spirals of various shapes were observed both on the ground and in 
microgravity for small particles of different sizes, shapes, and densities [9l[l0]. However, the 
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FIG. 2: Self-assembling particulate structures, (a), Azimuthal (left) and axial (right) snapshots 
of the spiral structures observed in the experiments of D. Schwabe et al. [lOj. The illuminated 
particles are seen as white specks. The closed spiral rotates in the direction denoted by the arrow, 
(b), Corresponding views of the structures obtained in our direct numerical simulations [22j. The 
particle instantaneous positions are denoted by dots. On the left plot the particles positions are 
superimposed on the plot of the instantaneous temperature field taken at the cylinder mid-plane. 
Here m — 2. 

coherent structures were observed only for certain parameter ranges. 

As the convective time-based Stokes number St = 2/9 (a/L)'^ Re ^ 10~^ ^ 10~^ <C 1, 
the particle dynamics is dominated by the viscous drag force. Here a is the particle size, 
L is the characteristic flow length scale, Re = UL/u is the Reynolds number, U is the 
characteristic flow speed, and u is the fluid kinematic viscosity. In previous experimental 
studies of the particle self-ordering researchers looked for additional forces and accounted 
for flow features speciflc to thermocapillary and free-surface flows [ilOj. However, we claim 
that this phenomenon is generic for the class of volume-preserving flows deflned above, and 
that an interplay of the particle inertia and the viscous drag force alone may cause the 
ordering. Therefore we anticipate this effect in such periodic flows, which abound in nature. 
Their further examples include flows in laminar stirred tanks p^, pyroclastic surges, and 
microfluidic flows ^18j. The arising singular spatial distributions of particles may have 
profound consequences for particle aggregation and transport; understanding and control of 
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the arising structures may be important as a technological tool. 

The starting point of our analysis is the reduced version of the Maxey-Riley equation of 
particle motion [19] , which can be written as 



Here t is time, v is the particle velocity, u is the fluid velocity, p is the ratio of the particle to 
the fluid density, r is the particle relaxation time, and g is gravity. The material derivative 

is taken moving with the fluid velocity at the current location of the particle. By the 
Stokes law, r = (2/9)a^/z/. This equation describes the Newtonian dynamics of small 
particles dominated by the inertia (including the added mass effect) and viscous drag forces. 
It neglects a reverse influence of particles on the flow and particle-particle interactions and 
is valid when a/L <C 1 and Ua/u <C 1. (While in an actual physical system particles may 
experience other types of hydrodynamic forces, our goal is to demonstrate that even the 
current 'minimal' model can produce the effect.) 

We performed numerical simulations, in which the particle dynamics governed by ([T]) was 
coupled to the full system of the Navier-Stokes equations describing the thermocapillary fluid 
flow [20l El]. We flnd, in physically realistic regimes, that particles assemble in dynamic 
spirals that closely resemble the experimental results (Fig. [2](b) and Supplementary Movie 
2 [i7j). Besides, we observed formation of particulate spirals that have the number of turns 
/ different from the wave mode m (also integer due to the cylindrical geometry) [22j . Similar 
to experimental flndings, each coherent structure is robust in a limited range of governing 
parameters. As soon as a parameter leaves the range, the coherent structure will disperse. 
We will come back to these features below. The major lesson learned from the simulations 
is the basic fact that a mere interplay of the inertia and viscous forces acting on individual 
particles can lead to the particle ordering. 

At this point one could suppose that the effect is speciflc to the thermocapillary flow. In 
order to demonstrate that that this is not the case, we study particle ordering in a generic 
analytical model of a flow with a rotating wave. For a cylinder of unit radius we assume: 
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FIG. 3: Distribution of particles in the equatorial plane rotating with the wave after n iterations 
of the Poincare map. (a) n=0, (b) n=150, (c) n=300, (d) n=3000. Here m = 3, i7 = 1, 6 = 0.3, 
= 0.5, and r{p — 1) = 10~^ (see the second term in Eq. [s]). In the physical space this picture 
corresponds to formation of a two-dimensional toroidal coherent structure that later transforms 
into a one-dimensional closed spiral. The projection of the spiral is showed in panel (d) by the red 
line. The physical formation time is about 10^. 

u^^\(f)^r^z) = — (3r^ — 4r^) sin (m(/)) sin 



u, 



H 

(1)/ , X / 9 ctx , . TTZ 



,r, z) = (r^ — r^) mcos (m(/)) sin — . (2c) 
Here H is the cylinder height, m is the (integer) wavenumber, and h is the wave amplitude. 



The equations (2b) mimic the Poiseuille flow in a flnite cylinder, and the equations (2c) - a 
rotating wave. The flow is scaled by [/, so that the maximum value of is 1 at z = H/2 and 
r = 0, or 1. This model is not intended to be a rigorous description of the flow observed in 
experiments or numerical simulations; rather, it is a phenomenological model that captures 
the essential features of flows with rotating waves. It satisfles impermeability boundary 
conditions. A model satisfying no-slip conditions on the solid walls and zero axial vorticity 
component on the free surface would be signiflcantly more complicated. Nevertherless, this 
simplifled model successfully reproduces the principle features of the particle ordering, as is 
showed below. 

Instead of solving ([!]), we deal with the 'inertial equation' obtained as the flrst order 
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FIG. 4: Synchronization of a particle and the wave. Here m = 3. (a) Schematic trajectory of an 
individual particle in the laboratory frame. (b),(c) The equatorial plane in two instances of time 
separated by the particle turn-over time Tp. The phase locking means the particle phase relative to 
the wave is the same for the consecutive locations '1' and '2'. If the particle gets out of phase, '2^', 
interaction with the flow structure will adjust the particle azimuthal drift (red arrow) to resume 
the phase. 

approximation of ([T]) in iSt: 



This equation describes particle motion as advection plus a perturbation due to the inertial 
effects. Since V -i^ = 0, it is the nonzero divergence of the second term that must be respon- 
sible for the phase volume changes accompanying formation of the accumulation structures. 



from six to three dimensions (in the frame rotating with the wave). Then from the viewpoint 
of dynamical systems the coherent structures are attractors in the three-dimensional space 
and can be readily studied by means of the Poincare section. We define the latter as the 
instances when particles cross the equatorial plane z = 1/2 (in the reference frame rotating 
with the wave). 

The model ([2]), ^ contains five parameters 77, m, a, and (p — l)r, and can be studied 
numerically. At this point our goal is two- fold: first, to show that the model may account for 
the spontaneous formation of one-dimensional structures from the initially three-dimensional 
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(3) 



The chief advantage of using ^ instead of ([T]) lies in the reduction of the particle phase space 



particle distributions for a realistic choice of parameters. Second, to demonstrate that it 
requires no fine-tuning of the parameters, i.e. the structures are stable for certain parameter 
ranges. 

We choose i7 = l,m = 3, and b = 0.3 to account for the flow and the rotating wave 
characteristics observed in experiments and our numerical simulations, and (p — l)r = 10~^ 
to account for the particle properties. The main tunable parameter is Q and for now we set 
it to 0.5. Numerical solution of ([2]), ([3]) bears out formation of coherent structures similar to 
the ones observed in experiments P [lO] and direct numerical simulations [21j. We clearly 
observe that particles do not need to touch the fluid-gas interface in order to align. Fig. |3] 
demonstrates that the coherent structures are stable flxed points of the Poincare map. It also 
shows that the particle ordering proceeds via two distinct steps: flrst, particles concentrate 
in the center and near the circumference of the Poincare section. Second, the cylindrical 
symmetry is broken and the particles are attracted to the stable flxed points. In the physical 
space these steps correspond to clustering of particles in two-dimensional toroidal coherent 
structures, and to transformation of the latter into one-dimensional closed spirals. The two 
processes have different characteristic times but take place simultaneously [22]. (These steps 
of the particle self-ordering had been noted in experiments [23j.) As a result, the motion 
of an individual particle in the spiral gets synchronized with the wave, and the resulting 
structure appears to rotate along with the wave in the laboratory frame. 

We flnd an explanation for this behavior in the nonlinear phenomenon of phase locking. 
Phase locking, alternatively called frequency locking or entrainment, typically occurs in 
dissipative systems with a weak interaction between modes with close frequencies [24j. In 
our case the turn-over particle motion is synchronized with the rotating wave oscillations. 
The physical essence of the mechanism lies in the adjustment of the azimuthal particle 
displacement after every particle turn-over 0^+1 — (j)^ due to the inertial interaction with the 
wave, see Fig. [4[ Consequently, the particle azimuthal drift Up^cfy is also modifled. 

It is a signature of phase locking that the ratio of frequencies, called the winding number, 
remains a constant rational over a range of control parameters [24]. In our case the winding 
number W = Cl/up^ where Q is the frequency of the wave oscillations experienced by the 
moving particle, and ujp is the frequency of the particle turn-over motion. Fig. [5] shows 
that while ujp changes as Q is varied over the range [0.39, 0.67], W remains identically equal 
to 1. The values of W change abruptly outside the phase locked regime and the coherent 
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FIG. 5: The winding number plateau. The winding number W (sohd, black) and the particle 
turnover frequency Up (dashed, blue) vs. the wave angular velocity fi. While ujp changes, W 
remains constant for a range of values of ft. Constancy of winding number is a signature of phase 
locking [24j . 

structures loose stability. We find other rational values for the winding number W = p : q 
for different parameter values and in certain regimes in direct numerical simulations [22j. 
The latter equality can be re- written as the 'resonant condition' pujp = qCl. (A similar 
'resonant condition' was noted in [lOj. Observe, however, that unlike resonance, phase 
locking does not require fine-tuning but provides it.) We note that Q = m{Q—ujp^cfy) ('angular 
Doppler shift'). Thus, the phase locking adjusts the particles azimuthal drift to maintain the 
'resonant condition.' The three integers p, g, and m fully determine the geometrical shape 
of the particulate spirals [22j . 

Synchronization due to phase locking is ubiquitous in nature. The present modeling 
suggests that PAS formation in thermocapillary fiows is another instance of this general 
phenomenon and that formation of one-dimensional coherent particulate structures should 
be encountered in other oscillatory vortical fiows when 1) the particle turnover motion 
is transversal to the direction of wave propagation and 2) the frequencies of the particle 
turnover motion are commensurate with the oscillation frequencies. 
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